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Abstract 

We perform a systematical analysis of the spin Hall effect (SHE) in the Cu alloys doped with a series of 5d elements, 
by the combined approach of density functional theory and Hartree-Fock approximation. We find that not only the 
spin orbit interactions (SOI) in both the 5d and 6p orbitals, but also the local correlations in the 5d orbitals of the 
impurities, are decisive on the sign of the spin Hall angle (SHA). Including all of these three factors properly, we 
predict the SHA for each alloy in the series. The signs of Cuir and CuPt are sensitive to perturbation of the local 
correlations. This observation is favorable for controlling the sign of the transverse spin Hall voltage. 
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1. Introduction 

Spintronics is a rich research field, not only for 
the wide applications of low energy consumption and 
energy transformation devices, but also for the new 
physics with the interplays among charge, spin, orbital, 
heat and so on 111]. The spin Hall effect (SHE), which 
converts the injected longitudinal charge current into 
the transverse spin current via the spin-orbit interaction 
(SOI), is crucial for the development of spintronic de¬ 
vices. The SHE is characterized by the ratio between 
the transverse and the longitudinal resistivities (or con¬ 
ductivities), called spin Hall angle (SHA). The magni¬ 
tude of SHA describes conversion efficiency between 
the charge current and the spin current, while the sign 
distinguishes the scattering direction of electrons, i.e., 
clockwise or anticlockwise into the transverse direction. 
It is well known that to raise the magnitude of SHA is 
the key to enhance the whole efficiency of the devices 
based on the SHE However, the sign of SHA 

is not yet utilized as a degree of freedom in spintronic 
devices. 

In the experiment of the dilute CuIr alloys, the dom¬ 
inant contribution to the SHE was verified to be an ex¬ 
trinsic skew scattering mechanism, and the SHA was 
measured to be positive 2.1% 0). We have found out in 
theory that including the local correlation effects of the 
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5d orbitals of Ir is decisive to obtain a positive sign of 
SHA consistent with experiment, and a small change of 
the 5d electron number of Ir with the local correlations 
may change the sign Si. 

On the other hand, among the previous theoretical 
works about the SHE of the Cu alloys doped with 5d el¬ 
ements, Fert and Levy estimated the SHA by the atomic 
electron numbers of 5d orbitals with SOI, obtaining a 
sign change of SHA in the middle of the series of 5d 
impurities li. In contrary, based on the ab initio cal¬ 
culations, including the SOI in both 5d and 6p orbitals, 
Fedorov et al. obtained an uniform sign of SHA con¬ 
tributed by the skew scattering among the series of 5d 
impurities from Lu to Pt Isll . 

In the present work, we analyze the SHE among the 
dilute Cu alloys with a series of 5d elements as impu¬ 
rities, and find the key factors which are decisive on 
the sign of SHA. Furthermore, we screen out the alloys 
whose signs of SHA are sensitive to perturbation of the 
local correlations as Cuir, which will be favorable for 
the sign control of SHE. 

2. Theoretical Approach 

The extrinsic SHE is caused by the spin-orbit interac¬ 
tions (SOI) in the orbitals of the impurities. As one of 
the mechanisms of extrinsic SHE, the skew scattering 
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generates the spin Hall resistivity linear with the impu¬ 
rity concentration in the dilute alloys. Based on the An¬ 
derson model 10], the contribution of skew scattering on 
the extrinsic SHE in the nonmagnetic Cu alloys with di¬ 
lute 5d impurities can be described by a single-impurity 
multi-orbital model, including the SOI in both the 6p 
orbitals ^ and the 5d orbitals ^ of the impurity with the 
parameters Ap and Aj, respectively, together with the lo¬ 
cal correlations of the on-site Coulomb repulsion U (U') 
within (between) the 5d orbitals, and the Hund coupling 
J between the 5d orbitals of the impurity |@]: 
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In Eq.Q, Cak is the energy band a of the Cu host, 
is the energy level of the orbital jS of the 5d impurity, 
and V(a,ff(k) is the hybridization between the orbital p of 
the impurity and the band a of the host. In Eq.©, the 
local correlations are included only within the localized 
5d orbitals f of the impurity, employing the relations of 
U ^U' + 2J^ and 7/t/ = 0.3. 

In the orbital with the orbital angular momentum /, 
the SOI split the states of the orbital into two groups of 
degenerated states with the total angular momentum j - 
I + 5 and the degeneracy Di± -2] + \ .\n the Cu alloys 
with the 5d impurities including SOI, according to the 
the Eriedel sum rule dQ, the phase shifts 6f can 
be calculated from the occupation numbers Ni± of the 
spin-orbit split states of the impurity and the Cu host. 


6f -Nf:)IDi^, (3) 


where l-l corresponds to the 6p states of the impurity 
and the 4p states of Cu, 1=2 corresponds to 5d states of 
the impurity and the 3d states of Cu, respectively. The 
total occupation number of the 5d states of the impurity 
can be obtained by 

( 4 ) 


The occupation number of each of the degenerate 5d 
states of the impurity nd± can be obtained by 

nd± = N'Z‘’/D2^, (5) 

where will be between 0 and 1. The occupation 
numbers of the impurities are defined via projections of 
the occupied states onto the Wannier states centered at 
the impurities as point defects and extended in the whole 
supercell. In the dilute alloys, a point defect cannot be 
charged in metal. Thus the total occupation numbers of 
the 6s, 6p and 5d states of the impurities are conserved 
to the atomic valence electron numbers Si, 

Nr+ n;'p+N fp = ( 6 ) 

where the constant C"”^=4, 5, 6, 7, 8, 9, 10, 11 and 12, 
for the impurities of Hf, Ta, W, Re, Os, Ir, Pt, Au and 
Hg, respectively. 

Based on the model in Eq.([T]i, by the Hartree-Fock 
(HF) approximation, the 5d states of the impurities can 
be considered as virtual bound states with a width pa¬ 
rameter A, and the 5d± states correspond to the ener^ 
levels of eo,d±, respectively. It has the relations of Jalst] 


Acot(OT!rf±) = eo,rf± - = £'o,d±, (7) 


where the er is the Fermi level and £’o,^/± are the energy 
levels of the 5d± states relative to the Fermi level. While 
the local correlations are included in the 5d states of the 
impurity as in Eq.(|2]i, it has the self-consistent relations 
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which will raise the 5d energy level and lower the 5d 
occupation numbers of the impurity SSi. The de¬ 
crease of N'r will induce the increase of Np'^’ accord¬ 
ing to Eq.®, which can be approximately estimated by 
fixing the ratios of I and The cor¬ 

responding phase shifts can be obtained by Eq.®. 

For the SHE contributed by the skew scattering in Cu 
alloys with 5d impurities, the SHA defined in terms of 
resistivity p llj] can be calculated from the phase 
shifts df of the p± channels and <5^ of the d± channels. 
Including the SOI in both the p and d channels, the SHA 
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© is obtained by |@] 
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(9) 

When the SOI in the p (or d) channels are omitted, with 
(5i = <57 — 6 ^ (or 62 — ^2 — 57 ), the simplified formula 
for SHA 0rf(bi,bi,(57,57) (or Qp( 6 ^, 6 ^, 62 , 62 )) can be 
obtained from Eq.®- 

Based on the model in Eq.© and the formula in 
Eq.®, to calculate the SHA for the Cu alloys with a se¬ 
ries of 5d impurities, we employ the combined method 
of the density functional theory (DET) and HE approxi¬ 
mation, as we have done before for the Cuir alloys |@] . 


3. Results and Discussions 


Eor the dilute alloys of Cu with a series of 5d impu¬ 
rities, including the SOI in both the 5d and 6 p orbitals, 
the 5d occupation numbers Nd and the SHA are calcu¬ 
lated by the DET-HHE method following Eqs.® and ®, 
without the local correlations (U-Q) or with the local 
correlations of U =0.5 eV in the 5d orbitals of the impu¬ 
rities, respectively, as shown in EiglT] Estimates from 
experiments give the parameter of the on-site Coulomb 
repulsion U of 5d orbitals in the range of 0-1.0 eV for 
pure metals lll4ll . 0-1.5 eV for insulating compounds 
1 1511 . less than 0.4 eV for the ferromagnetic metal of 
osmates |15|, about 0.5 eV for iridates ifl^ . and 0.5- 
0.6 eV for Pt alloys 114]. Below we take U-0.5 eV to 
compare the whole series of 5d impurities in Cu alloys. 

As U switches from 0 to 0.5 eV, according to Eq.® 
based on the Anderson model mss, the 5state will 
be shifted up, accompanied by a decrease of Nd- The 
amount of the shift varies among the 5d impurities, as 
intuitively shown in Eigl^ The Nd lowered by U for 
each of the 5d impurities are marked in EiglTfa). In the 
middle of the series with Nd around the half-filling value 
of 5, the Nd is obviously affected by U, corresponding 
to the schematic picture in Eigl^that the peaks of the 
5d band from Hf to Pt are relatively close to the Eermi 
level. 

Previous theories leave the following discrepancy re¬ 
garding the sign of the SHA in the Cu alloys with 5d 
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Figure 1: (a) The occupation numbers of the 5d orbitals of the im¬ 
purities Nd and (b) the SHA (©(d|, dj, d^)) of the Cu alloys with 
5d impurities under the local coiTelation U=0 (black triangles) and 
U=0.5 eV (red circles) in 5d orbitals, respectively. The experimental 
data of SHA in CuIr 0 is mai'ked by the red cross in (b). 
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Figure 2: Schematic picture of the total density of states (DOS) of the 
5d virtual bound states, for the impurities in the dilute alloys CuHf 
(blue), CuPt (red) and CuHg (black), with the on-site Coulomb re¬ 
pulsion U=0 and U=Q.5 eV, respectively. Sf is the Fermi level. The 
amount of the shift by U depends on the atomic number of the impu¬ 
rities. 
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Figure 3: (a) The occupation numbers of the 5d orbitals Nj of Pt (black curve), Au (red curve) and Hg (blue curve) in the Cu as functions of the 
local coiTelation U in 5d orbitals. The SHA of the (b) CuPt, (c) CuAu and (d) CuHg alloys as functions of U. In (b), (c) and (d), the functions of 
SHA 0(d|,(5^,(5j,(5p, &p( 6 ^, 6 ^, 62 ,Si) and are denoted by black solid, blue dash and red dash-dot curves, respectively. 


impurities contributed by the skew scattering. Fert and 
Levy gave a sign change of SHA in the middle of the 5d 
series, including the SOI only in the 5d orbitals without 
local correlations, and taking the phase shift of p chan¬ 
nel as a parameter 101 ■ Fedorov et al. obtained that the 
signs of SHA among a series of 5d impurities from Hf 
to Pt do not change, including the SOI in both the 5d 
and 6p orbitals but without local correlations. Their re¬ 
sult for Cuir is opposite in sign to the experiment under 
the consistent definition of SHA via resistivity |0]. 

In comparison, in FiglHb), the series of SHA under 
U-Q is very close to the ab initio results of the Cu alloys 
with the impurities from Ta to Pt in Ref.jst] (with the 
SHA consistently defined via resistivity), as expected. 
However, the series of SHA deviate obviously between 
U-Q and t/=0.5 eV, except for the CuHg case with the 
5d states nearly fully occupied. It is remarkable that for 
the cases of CuIr and CuPt, the SHA changes sign be¬ 
tween U-0 and U—0.5 eV. Moreover, with U-0.5 eV, 
Cuir gives a SHA of -1-1.6%, consistent to experimental 
value of -h 2.1% 101, as marked in Fig [Tib). These facts 


demonstrate that the SOI in both the 5d and 6p, together 
with the local correlations in the 5d orbitals of the impu¬ 
rities, are decisive on the sign of SHA of the Cu alloys. 
All of these three factors are necessary to be included 
properly to estimate the SHA contributed by the skew 
scattering. 

Fig|3] shows the changes of Nd and SHA along with 
U for the cases of dilute CuPt, CuAu and CuHg alloys. 
The results of CuPt in Figsj^a) and (b) are quite similar 
to the results of Cuir shown in Ref. 1^ . Close to the gen¬ 
erally realistic value of U-0.5 eV, both the Cuir 1@1 and 
CuPt have a sign change of SHA between U-Q.3 and 
0.4 eV, corresponding to the change of Nj from U =0.5 
eV by as few as about 0.3 electrons. In the experiment 
of the heavily doped CuPt alloys with 10% of Pt, the 
charge transfer into the 5d orbitals of Pt is 0.25 elec¬ 
trons IB. If it were possible to devise a way to change 
the Nd in Cuir or CuPt by about 0.3 electrons, or if the U 
could be altered by about 0.1 eV, the sign of SHE could 
be externally controlled. 

Within the range of U from 0 to 1.5 eV, the sign of 
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SHA of the Cu alloys with 5d impurities from Hf to Os 
is stably negative. There are more sign changes of SHA 
around U-1.25 eV for CuPt and CuAu, as shown in 
FigsI3b) and (c), but such a large U is unlikely to be rel¬ 
evant for 5d impurities. In contrast to CuPt, the changes 
of Nd and SHA in CuAu and CuHg are extremely slow 
with U, as shown in Figs[3la), (c) and (d). This is be¬ 
cause the 5d states of Au and Hg are nearly fully oc¬ 
cupied, and the center of 5d band is so deep below the 
Fermi level that the shift by U induces little change of 
Nd, as shown in Fig|2 Increasing U from 0 to 0.5 eV, 
the decrease of Nd of Au and Hg is as small as 0.13 and 
0.03 electrons, respectively. Also due to the nearly full 
occupation of the 5d states, the phase shifts of d channel 
are very close to zero and the spin-orbit splits between 
d+ and d- states are very small for CuAu and CuHg. 
As shown in Figsjjc) and (d), with the U from zero 
to around 0.5 eV, the SHA function 0(d|, 5^, is 

closer to ©p(d|, dj", 62 , 62 ), rather than Qd{d\, 61 , 62 , 5p, 
indicating that the skew scattering by the SOI in 6p or¬ 
bitals overwhelms the contribution by the SOI in 5d or¬ 
bitals. 

In summary, by the combined method of DFT-hHF, 
we analyze the SHE in the Cu alloys with a series of 
5d impurities. We stress that it is necessary to include 
the SOI in both the 6 p and 5d orbitals, together with the 
local correlations in the 5d orbitals of the impurity, to 
estimate the SHA of these alloys. The results reveal that 
the sign of SHA of Cuir and CuPt would be sensitive to 
perturbation of the local correlations, which is favorable 
for the sign control of SHE performing as a spin current 
switch. 
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